Families of unitary matrices 
achieving full diversity 



o 
o 

!^ 
< 



C/3 



> 

o 

00 

o 
in 
o 

O 



X 



Frederique Oggier 
Ecole Polytechnique Federale de Lausanne 
Laboratoire de mathematique algorithmique 
1015-Lausanne, Switzerland 
Email : frederique. oggier @ epfl. ch 



Abstract — This paper presents an algebraic construction of 
families of unitary matrices that achieve full diversity. They are 
obtained as subsets of cyclic division algebras. 

L Problem statement 

In the context of noncoherent multiple antennas channel 
coding, research has been done on constructing families of 
unitary matrices with full diversity, that is, satisfying that the 
determinant of the difference of any two matrices in the family 
is nonzero. Among the algebraic approaches to this problem, 
the theory of fixed-point-free groups and their representations 
has been exploited in [6], while representations of Lie groups 
has been investigated by Jing and Hassibi (see e.g. [1]). 

At the same time, division algebras for space-time coding 
have been introduced in the context of coherent MIMO sys- 
tems [5]. These algebras became of great interest, since they 
naturally provide a linear family of fully-diverse matrices. 

The aim of this work is to show that division algebras 
(in particular cyclic division algebras) can also be used to 
construct fully diverse unitary matrices. 

The paper is organized as follows. In the next section, we 
recall the basic facts about cyclic algebras. In section [111] 
we explain how the condition of being unitary for a matrix 
can be translated into first a constraint on an element of the 
algebra, and second a constraint on a commutative subfield 
of the algebra. This contains a constructive proof that yields 
a way of exhibiting unitary matrices. The whole process is 
illustrated in a worked out example in section Hvl 

Remark 1: In the following, we choose the dimension of 
the algebra to be 3 for the sake of simplicity. The same theory 
can be generalized for any dimension n. 

II. Cubic cyclic algebras 

In this section, we briefly recall what is a cubic cyclic 
algebra, and how it provides a linear family of 3 x 3 fully- 
diverse matrices. Let L, K be two number fields. 

A. The algebra structure 

Let L/K be a Galois extension of degree 3 such that its 
Galois group G = Gal{L/K) is cyclic, with generator a. 
Namely, G — {(J,a^,a'^ — Id}. Such an extension is called 
cyclic. Denote by K* (resp. L*) the set of non-zero elements 
of K (resp. L). We choose an element 7 e K*. We construct 
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a non-commutative algebra, denoted A = {L/K,a,-f), as 
follows: 

A = L®eL®e^L 

such that e satisfies 

= 7 and Ae = eo'(A) for X E L. 

Such an algebra is called a cubic cyclic algebra. It is a right 
vector space over L, and as such has dimension {A : L) — 3. 

Cubic cyclic algebras naturally provide linear families of 
matrices thanks to an isomorphism between the split algebra 
A ®K L and the algebra Mj,{L), the 3-dimensional matrices 
with coefficients in L. This isomorphism, denote it by h, is 
given explicitly. Since each a; G .4 is expressible as 

a; = xq + exi + 6^x2, Xi E L for all i, 

it is enough to give h{xi 1) and h{e (g) 1). We have that 

h:A(E)KL = M3{L) (1) 

is given by, for all i. 



XiCS)l 



Thus the matrix of ft, (a; ® 1) is easily checked to be 

xo ia-ix2) jo-'^ixi) 
xi a{xQ) ja'^{x2) 
X2 a-ixi) a'^ixo) 
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Remark 2: Notice that (|2jl is also the matrix of left multi- 
plication by X in the basis {1, e, e^}. 
We thus start with the family of matrices 

{/ Xq 7(7(2:2) jcr'^ixi) \ 
X= \ xi a{xo) -fcr'^ix2) \\xo,Xi,X2&L 
\ X2 a{xi) cr'^ixo) J 

which is linear, since it has an algebra structure. 



B. The diversity property 

Recall that the diversity product (^{C) of a set C of AI unitary 
3x3 matrices Xi, . . . , Xm is the minimal diversity distance 

C(C) :-imin|det(X,-X,)|i/3. 

A set of matrices with C(C) > is said to have full diversity. 
If C is linear, then the above definition simplifies to 

C(C) := l min |det(X)|i/3 

in which case full diversity is obtained if all matrices are 
invertible. Take now C C A, A an algebra. If furthermore 
^ is a division algebra, all matrices in C are by definition 
invertible. Thus C is fully diverse. Summarizing, if there is C 
a subset of unitary matrices in ^ a cyclic division algebra, 
then this family will automatically be fully diverse. 

To decide whether a cyclic algebra is a division algebra, the 
following criterion is available: 

Proposition 1: [3, p. 279] Let L/K be a cyclic extension 
of degree n with Galois group Gal{L/K) =< cr >. If 7 and 
its powers 7^, . . . , 7""^ are not a norm, then {L/K, a, 7) is 
a division algebra. 

III. The unitary constraint 

Suppose ^ is a division algebra. We now have a linear 
family of invertible matrices among which we are looking for 
unitary matrices, i.e., X G C such XX* = I3, where * denotes 
the transpose conjugate. 

A. The unitary constraint in the algebra 

We take advantage of the matrices coming from the algebra 
A and use the following correspondances: 

X ^ A'-^ X ®l e A®K L ^li. ^ MziL) 

We thus translate the condition of "being unitary" for the 
matrix X into a condition on the element x in the algebra. 
We will show that A can be endowed with an involution a, 
and that 

XX* =13 h{x®l)h(x®l)* = 1 ^ xa{x) = l. 

Remark 3: Defining an involution on the algebra (i.e., a 
map that satisfies the three properties described in the proof 
of Proposition |3 and checking that it is well defined on the 
algebra of matrices is technical, but this formalism is required 
to be sure that our objects are well defined. This subsection 
ends with an example that illustrates the theory. 
Let us first define an involution on A. 

Proposition 2: Let : L — > L be an involution on L 
such that aL commutes with all elements of Gal{L/K). Let 
z = 70;^ (7) and a : A A such that 

a{xo + exi + e'^ X2) = aLixo) + e~^za'^{aLixi)) 

+e^'^z'^a~'^{aL{x2)). 

Then a defines an involution on A. 



Remark 4: Note that the condition that ul commutes with 
all elements of Ga\{L/K) implies that aL{K) = K. Indeed, 

cr{oiL{k)) — aL{(j{k)) — aL{k), for any k £ K, 

showing that aL{k) is fixed by a. In particular, z & K. 
Proof: Note first that we have ea{e) — z. Check that 

1) a{x + y) = a{x) + a{y) for all a;, y e A 
This is clear 

2) a{e^ yjC^Xi) — a{e^Xi)a{e^yj) for Xi,yj G L. 
If i + j < 3, we have 

a{e^yje'xi) = a{e'+' (T\yj)xi) 

^ e- V- (a^ (a;, ) (a' (y, ) ) ) 

= e- z'+^a-^'+^\aL{x,))a-={aL{y,)). 

Now, the right handside term is given by 

a{e'xi)a(e^yj) = er'z'a^''{aL{xi))e^^ z^a^^ (aLiyj)) 

= e-(^+-'')z'+-'a-(*+^)K(x,))a-^(aL(yj)). 

Similarly, if i + j > 3, i + j = 3 + fc, < fc < 2 and 
the same computations hold. 

3) a{a{x)) = X for all x ^ A. 
We have 

a(a(e'x,,)) = a{e~^z'a~\aLix,))) 

= aL{(T~'{aL{x^))aL{z')a{e-') 
= a {xi)aL[z )z e . 

Since z is fixed by ul, we get that 

a{a{e^Xi)) = e*Xi. 

■ 

The involution a defined in the above proposition is extended 
to the spHt algebra A ®k L = M^iL) as follows. 

a® ul A ®K L ^ A ®K L. 

It is used to define an involution ah on Ai^lL) via the 
isomorphism h: 

ah — h o {a ® ul) o h'^ ■ (3) 

Proposition 3: Let X — h{x®l) and z be as in the hypoth- 
esis of the previous proposition. If z = 1, then a;i(X) = X*. 
Proof: Recall first that 

(xq j<7{x2) ja'^ixi) \ 

xi a{xo) 'y(T^{x2) 

X2 (Jixi) a^{xo) J 

We have 

ah{X) = ah{h{x®l)) 

= ho {a® aL){x ® 1) 

= h{a{x) ® aL{l)) 

= h{a(x) ®l)aL{l)- 



Recall that e^^ = 7^^e^, 7"^ = z'^aLij) ^"'1 that h{a{x)^ 
1) is the matrix of multiplication by a{x) (see Remark |2}. 
Since 

a{x) = aLixo) + e~^za^^{aLixi)) + e^'^z'^a^'^{aLix2)) 
= aL(xo) + e"f~'^z^a{aL{x2)) + e^^^^ za'^{aL{xi)), 

we have 

h{a{x) (g> 1) = 

aL{xo) zaL{xi) z'^aL{x2) \ 

ai (7) zcr(aL (2:2)) cr(aL(xo)) zcr(aL(xi)) . 

(7)0'^ ("1,(2:1)) aL(7)2:o-^(aL(a;2)) cr^(aL(a;o)) / 

Since z = 1, and a commute, and is multiplicative, 
we get the desired result. ■ 
Notice that clearly 2 = 1 is the only possible choice for the 
involution to be the conjugate transpose. These formal 
proofs finally yield the desired result: 

Corollary 1: We have the following equivalence: 

XX* = I3 xa{x) ~ 1. 

Example 1: Let K = Q(C3) be a cyclotomic field, where 
(■3 is a primitive third root of unity, and let L ~ KQ{9) be 
the compositum of K and a totally real cubic number field 
Q{0), with discriminant coprime to the discriminant of K and 
cyclic Galois group Gal((Q(6')/(Q) =< a > (see Figure[i]i. We 
consider the algebra A = (L/K, a, 7), where 7 = C3. 

The involution on L is given by 

aL ■ L L 

ao + ai9 + a29'^ ^ T{ao) + T{ai)e + T{a2)e'^ 

where t is the generator of the Galois group Gal(Q(C3)/Q). 
Namely, r(6o + 61(3) = ^'o + &iCl and = -(3 - 1. The 
involution satisfies that r commutes with a, so that the 
involution 

a: A ^ A 

x = xq + exi + e^X2 ^ a{x) 

where a[x) = aL{xo) + eC,la{aL{x2)) + e'^C,l(7'^{aL{xi)) is 
well defined by Proposition |2| 
We have 

007\/0 10\ /100\ 
1 1 = 1 

010/\700/ \001/ 

<;=4> 77 = 1 <;=^ ea{e) = 1. 

B. The unitary constraint in commutative subfields 

We now show how the problem of finding unitary elements 
in the algebra A can be reduced to find elements of norm 1 
in commutative subfields of A. 

Proposition 4: Let A = {L/K,a,j) be a cyclic division 
algebra, and x A such that x ^ ±1. The following 
statements are equivalent: 

1) xa{x) = 1. 

2) There exists u ^ A* such that ua{u) = a{u)u and 
X = ua{u)^^ — a{u)^^u. 

Proof: The sufficient condition is clear. Let us prove the 
necessary condition. Assume that xa{x) = 1. Let M denote 



A 



L 




Q 



Fig. 1. The cyclic algebra A= {L/KyCr,^). 

the subfield of A generated by K and x. It is commutative and 
satisfies that a{M) — M, since a{K) = K and a{x) ~ x^^. 
Thus 

= {y G M I a{y) = y}, 

the subfield of M fixed by a is well-defined. 

Claim: M/M°' is a quadratic extension with Galois group 

Gal(M/M") = {/c;M,a|M}. 

It is enough to prove that a\M 7^ IdM. If oi\m = Id-M, 
then xa{x) = = 1, implying that x = ±1, which is a 
contradiction. 

The condition xa{x) = 1 becomes Nj^j^M^i^) = 1. By a 
corollary of Hilbert 90 Theorem, there exists u £ M* such 
that X = u/ a{u). ■ 
The above proof gives a way of building unitary elements of 
the algebra A. Take a commutative subfield M of A such that 
a{M) = M but with y e M such that a{y) ^ y, so that M" 
is not M itself. Take u G M* and compute x = u/a{u). The 
element x & A will satisfy xa{x) = 1. 

The next step is thus how to build commutative subfields 
of A which are stable by a. 

Definition 1: [4][p. 113] Let ^ be a cubic cyclic algebra. 
For X ^ A, define its reduced characteristic polynomial Xx as 
the characteristic polynomial of h{x ® 1). 
Let X € A. Its reduced characteristic polynomial is given by 

(Xo-X -fCF{x2) l(j'^{xi) \ 

xi (j{xo)-X -ra^{x2) 
X2 crixi) cr^ixo) - X ) 

It is shown that for each x € A, its reduced characteristic 
polynomial lies in K[X] [4] [p. 113]. 

Lemma 1: If x ^ K, the center of A, then Xx is irreducible 
over K. 

Proof: Either Xx splits into three linear terms in 
which contradicts the hypothesis that x K, or it factors into 
one linear term in K[X] and one quadratic term irreducible 
over K. But {A : K) = 9, so that A cannot contain a quadratic 
subfield. ■ 
If X ^ K, K[X]/ {xx{X)) is a commutative subfield of A, of 
degree 3 over K. 
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-X 








1 




7- 







Example 2: Let A {L/K, a, 7) and ^ = L ® eL © e^L. 
Let Xe be the reduced characteristic polynomial of e, given by 



XeiX) = det 



By Lemma^] Xe is irreducible and K[X]/xe{X) = i^(e) is 
a commutative subfield of ^. 

Remark 5: We have not discussed whether the commutative 
subfields that we build contain a quadratic subfield fixed by 
a. This will be illustrated later in the example of Section Hvl 

IV. Families of fully-diverse unitary matrices 

In this section, we consider a particular cyclic division 
algebra and show how to use it to build families of fully- 
diverse unitary matrices. 

Let K = QiCs) and L = Q(C7+Cf \ Ca) be the compositum 
of ^(Cs) ™d Q(C7 + C'T^)' '■^^ maximal real subfield of the 
cyclotomic field <Q(C7) (see Figure Q with — (^^ + (7^)- 
We have Gal(L/Q(C3)) = (ff), with a : Cr + Cr^ ^ Cr+Cr^- 

Let A = (Q(C7 + C7\C3)/Q{C3),<J, Cs) be the correspond- 
ing cyclic algebra. This is a division algebra [2]. As already 
explained in Example [l] the involution a on ^ is given by 

a : A A 

X = Xq + exi + e'^X2 I— > a{x), 

aL{xo) + eC^aiaLix2)) 

/-2 

^3 

Note that is here given by the usual complex conjugation. 

A. Commutative subfields of A 

We consider the construction of commutative subfields of 
A which are stable by a. The first obvious subfield of A 
one can think of is L. By definition, it contains a totally real 
quadratic subfield, namely Q(C7 + C7^)- Then, as explained in 
Example 121 we can consider K{e), with minimal polynomial 
Xe{X) ^ X^ - C3. Thus K{e) = QiCg), with maximal real 
subfield Q(C9 + Cg^) (See Figure|2}. 

Let us now try to determine more systematically which are 
the commutative subfields AI of A such that M" ^ M. 
Clearly Af" <Z {x ^ A \ a{x) = x}. We thus look for 
conditions so as to satisfy x = a{x). 

Lemma 2: Let x — Xq ^ exi + 6^X2, with Xi G L, that is 
Xi ^ Vt + CsWi, Vi, Wi G Q(C7 -I- Cf ^) for i = 0, 1, 2. We have 

{Xq = aLixQ) 
Vl = -<t(u2) 
Wl ~ a{w2) + vi 
Proof: This is a straightforward computation. Identify 
the coefficients of the power of e 

aLicr{x2))^^^ 
aL(CT2(a;-l))7-i 



where a{x) 
aLibo + hCs) = bo + biC^ 



e^(^^a^{aL{xi)) and 



then develop and using that 7 = ^3, 
and the coefficient of (3. 



identify the constant term 





minimal polynomial of u 


discriminant 


u = e + {l + (3)e-e^ 
j/ = 29+(l + C3)e-e2 
v = 3e+(l + C3)e-e2 
i/ = 4e+(l + C3)e-e2 
v = 5e+(l + C3)e-e2 


+ - 5X - 3 
X3 - X2 - 12X + 1 

-6X -I 

~ IIX + 9 
X3 - X2 - 61X - 13 


2^ ■ 3 ■ 47 
11 • 659 
33 • 31 
3137 

2^ ■ 307 ■ 727 



TABLE I 

Examples of commutative subfields of A. 





Fig. 2. The algebra A and some of its commutative subfields 



Example 3: Let x = xq + exi + 6^X2 S A, with Xi = 
Vi + C3^j' ''^i' Wi € Q for i ~ 0,1,2. The conditions of 
Lemma|2]are vi — -~V2 and wi = W2 + vi. This defines the 
number field Q{(q + Cg^^), the maximal real subfield of the 
cyclotomic field ^(Cg)- Indeed, we have that y e Q(C9 + Ccr^) 
can be written as follows 



2/= 2/o + yi(C9 + C9"') + 2/2(C9 + C9"')', y^eQ, i = 0,i,2 
= yo + yi(C9 ' C9 - Cl) + y2(2 - C9 + C9 - Cg) 
= {yo + 2^2) + [(yi - y2) ~ y2C3]C9 + [(2/2 - yi) - yiC3]C9 

This finds "formally" a commutative subfield that we already 
found "naturally". 

Other examples can be found in Table |l| 



B. Unitary matrices in A 

We now illustrate how to build unitary matrices in the 
commutative subfield Q(Cg)/Q(C3) of A that we built in the 
previous subsection. 

Take for example the element 

X - l + Cg + Cg'+Cg eQ(C9) 
= (l + C3) + e + e2C3 

As a matrix, x can be represented as 



X = 



1 + C3 CI 

1 

C3 




We have 

a{x) = -Cg'-Cl + Cg'-Cl eQ(C9) 

Again, as a matrix, a{x) can be represented as 

/-Ca 1 CI \ 

Cs -C3 1 

V CI C3 -C3 ; 

which can be checked to be X*. We have 

x/aix) = 1/19(-10 + I6C9 + Cl - 4C| + 14C| + 8C|) 

which has norm 1, so that by Corollary [0 the matrix X(X*)^^ 
is unitary. This can be easily verified, since 

X(X*)-i ^ 

-0.421 - 0.182i 0.473 + 0.638i -0.157 + 0.36i V 

-0.236 - 0.319Z -0.421 - 0.182i 0.473 + 0.638i 

-0.789 + 0.09i -0.236- 0.319i -0.421 - 0.182i / 

Notice that this procedure can be applied 

• to any element of Q(C9)/Q(C3)' ™d for each, it will 
give a unitary matrix. There may obviously be some 
redundancy. Typically, if the element x is invariant by 
a, the above procedure will yield the identity matrix. 

• to any commutative subfield of A, assuming that it 
has a quadratic subfield fixed by the involution, so that 
we obtain a family of unitary matrices for every such 
commutative subfields of the algebra. 

• to any cyclic division algebra one starts with. 



V. Conclusion and future work 

In this paper, we showed how to build fully diverse matrices 
using cyclic division algebras. The main idea was to translate 
the condition of being unitary for a matrix into another 
condition that applies on commutative subfields of the algebra. 
We also showed that given a cyclic division algebra, several 
families of unitary matrices can be computed. Unlike the 
property of being unitary, the full diversity naturally came by 
choosing A a division algebra. 

In this paper, we focussed on the existence of unitary 
matrices in division algebras. We did not investigate further 
the properties of these matrices. In particular, we have not 
discussed yet their applications to noncoherent channel coding. 
Namely, we have not given explicit codebooks and analysis of 
their diversity. However, recent research work in that direction 
is already yielding promising results. In that perspective, there 
are several research topics we could suggest, among which: 
general lower bounds on the diversity, classification of the 
codebooks one can obtain from these families of unitary matri- 
ces, and in particular, criteria to design the "best" codebooks. 
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